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ABSTRACT 

Whole atmosphere seeing /3q is the most important parameter in site testing measure- 
ments. Estimation of the seeing from a variance of differential image motion is always 
biased by a non-zero DIMM exposure, which results in a wind smoothing. In the pa- 
per, the wind effects are studied within short exposure approximation, i.e. when the 
wind shifts turbulence during exposure by distance lesser than device aperture. The 
method of correction for this effect on the base of image motion correlation between 
adjacent frames is proposed. It is shown that the correlation can be used for estima- 
tion of the mean wind speed V2 and atmospheric coherence time to- Total power of 
longitudinal and transverse image motion is suggested for elimination of dependence 
on the wind direction. Obtained theoretical results were tested on the data obtained 
on Mount Shatdjatmaz in 2007-2010 with MASS/DIMM device and good agreement 
was found. 
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1 INTRODUCTION 

The measurement of differential image motion by using Dif- 
ferential Image Motion Monitor (DIMM) is widely spread 
method for characterization of the atmospheric optical tur- 
bulence (OT) (Sarazin & Roddier 1990; Vernin & Munoz- 
Tunon 1995). The differential image motion corresponds to 
the fluctuations of the difference of wavefront tilts (or differ- 
ence of angles-of-arrival) in two close apertures. The DIMM 
is simple and robust instrument which suits well for long- 
term and field campaigns. Interpretation of its output data 
based on the theory of light propagation through turbulent 
media is quite straightforward: variance of the differential 
image motion is unambiguously connected with seeing fio 
(Sarazin & Roddier 1990; Martin 1987). 

However, there are several effects introducing biases and 
random errors in DIMM results (Tokovinin 2002; Tokovinin 
& Kornilov 2007). One of them arises due to the fact that 
CCD camera captures star images which are used for cen- 
troids estimation with non-zero exposure. During typical ex- 
posure time r ~ 10 ms, a wind shifts the turbulence by 
distance comparable or more than aperture D. Therefore 
estimated centroids inevitably undergo the wind smoothing. 

Necessity of reduction for this effect led to appearance 
of both theoretical (Martin 1987) and experimental (Soules 
et al. 1996) studies of impact of the exposure on power of 
the differential image motion measured with DIMM. The 
methods of correction developed in these works require an a 
priori information about wind or some temporal characteris- 
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tics of the image motion estimated during the measurement 
process. Usually DIMM data is corrected for this bias using 
interlaced exposures. This method consists of alternation of 
single and double exposures (for instance 5 and 10 ms). See- 
ings are measured separately for these two sets and final es- 
timation calculated as certain combination of them (Sarazin 
1997; Tokovinin 2002). 

Realization of the measurements with modern CCD 
cameras changed situation significantly. It became possible 
to use shorter exposures and faster frame rates. On the other 
hand, the method of the interlaced exposures can be hardly 
used with these cameras. 

Alternative solution of the problem consists in a mea- 
surement of covariance of the image motion between adja- 
cent frames, which is defined by characteristic timescales of 
the process. In this paper we describe this method and prove 
theoretically its applicability. Also, it is shown that the co- 
variance can be used to estimate some effective wind in the 
atmosphere. 

In two first sections, we consider general expressions 
defining power of differential image motion and effect of 
wind smoothing for a single turbulent layer. The dependen- 
cies of the power (variance of differential image motion) on 
wind speed and wind direction and instrumental parameters 
are studied as well. 

In the section 4, we show that in short exposure regime, 
when the wind shear vt <C D, where v is the wind speed, 
the wind effects can be described as a quadratic functions of 
the shear. This approximation allows us to establish simple 
relation between the wind effects on variance and covariance. 

The section 5 generalizes obtained results for the case 
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of the whole atmosphere. Expressions of the wind effects on 
variance and covariance depending on some effective wind 
are given. Then we describe the correction of measured vari- 
ances to zero exposure and define an applicability and ac- 
curacy the method. Obtained formulae permit to evaluate 
square averaged wind Vi weighted with turbulence inten- 
sity and calculate the atmospheric coherence time To- This 
is demonstrated in the Section 6.3 on the base of analysis of 
real data. 

The last section contains discussion of the results and 
recommendations on optimal measurements and process- 
ing. Details of evaluation of needed integrals over two- 
dimensional spatial frequency are given in appendixes. 

The numerical estimates are computed for DIMM de- 
vice which we use for the OT measurement on Mount Shat- 
djatmaz (Kornilov et al. 20f0). Main parameters of the in- 
strument are aperture diameter D — 0.09 m, dimensionless 
baseline b — 2.18, exposure r = 4 ms before 2009 December 
13 and 2.5 ms after, period between exposures T = 5 ms. 



2 BASIC RELATIONS FOR SINGLE 
TURBULENT LAYER 

2.1 Spatial spectrum of the differential image 
motion 

Equations for spatial spectrum of the differential image mo- 
tion were given in classical works by Fried (1965) and Mar- 
tin (1987). In the paper of Martin (1987), these equations 
were developed by successive filtering of initial spatial spec- 
trum of light wave phase distortions on the assumption of 
Kolmogorov OT model. It is known that spectral power den- 
sity $ of phase fluctuations is proportional to the intensity 
A J = C 2 A/i in homogeneous and isotropic turbulent layer: 

'2tt\2. 
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where / is the modulus of the 2D spatial frequency, A is 
the wavelength and ro is the Fried parameter. For practi- 
cal purposes A J is more appropriate because this value is 
layer-by-layer additive and used for characterization of ver- 
tical OT profile. The phase spectrum is axisymmetric, i.e. 
depends only on the / and has the dimension m~ 2 . Eq. (1) 
is valid on condition that propagation effects of disturbed 
wavefront can be neglected (near-field approximation). 

Hereinafter we use coordinate system different from one 
used by Martin (1987). The x axis is placed along line con- 
necting centres of apertures. In frequency domain, coordi- 
nates of a point are defined by modulus / and position angle 
<f> being reckoned from the x axis. As usual two components 
of the differential image motion are considered: longitudi- 
nal is measured along the x axis (hereafter l-motion) and 
transverse is measured along y axis (t-motion). 

The spectrum F of the image motion can be represented 
as a product of phase spectrum, gradient filter, aperture fil- 
ter and differential filter (Martin 1987; Tokovinin 2002). For- 
mulae for these filters were transformed to used coordinate 
system. Gradient filter G for I- and t-motions has the form 
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where factor (A/27r) 2 is needed to pass from the slope in m -1 
to the angle- of- arrival in radians, measured in the method. 

Expression for the aperture filter A depends on how 
image centroids are evaluated. This issue is quite impor- 
tant because the wavefront cannot be considered flat within 
aperture of the diameter D and corresponding image in 
focal plane is neither diflractive nor even axisymmetrical. 
Normally two models are considered: gravity centre of the 
image (g-tilt, see Martin 1987) or position corresponding 
to the normal of plane approximating wavefront within 
aperture (z-tilt, see Fried 1975). The difference between 
these approaches was analyzed in the paper of Tokovinin 
(2002). Both filters are axisymmetrical and can be expressed 
through Bessel functions J n as 
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We will use A 9 (Df) by default since an exact form of aper- 
ture filter is not critical for analysis. 

The differential filter corresponding to two identical 
apertures placed at the distance B from each other along 
x axis is given by the expression: 



H(f, (j>) = 4sin 2 (7rB/cos ( 



(4) 



The final spectral power density F(f, 0) of the differen- 
tial image motion measured with DIMM is defined by the 
relation: 



FiAf, <t>) = *(/) Gi,t(f, 0) A(Df) H(f, 



(5) 



It is worth noting that after substitution of (1) and (2), the 
spectral density will not depend on wavelength A of incoming 
light. 

2.2 Power of differential image motion 

Variance of the differential image motion af t can be com- 
puted as integral of the Fi ;t (f, </>) over all spatial frequencies. 
Because integration is being carried out in polar coordinates, 
factor / has to be added. Let us group the spectral filters 
in order to separate terms depending and not depending on 
polar angle <j> for consequent averaging by this angle. Also, 
it is convenient to use the dimensionless frequency q — fD 
and the dimensionless baseline b = B/D. In this case, the 
variance of the differential image motion is 
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where the factor ^i,t{q) is a result of averaging of all terms 
depending on angle <f> over this angle. In case of the Z-motion, 
it becomes 

— — / sm {nbq cos <j>) cos (f>d(f> 

^ J -IT 

= 1-J (27rbq) + J 2 (2nbq), (7) 

and for the t-motion: 
2 
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bq cos (j>) sin (j>&4> 

= 1 - Jo (2-Kbq) - J 2 {2nbq) . (8) 
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For derivation of ^{q), the relations from Append. A were 
used. In cases when it is not important, we will discard sub- 
scripts designating component of the image motion. 

Sum of the variances of both components ct 2 = of + 
of (hereafter total differential image motion: c-motion) is 
very convenient for theoretical analysis and interpretation 
of observational data. Evidently, its gradient filter G c (f) = 
(A/) 2 is axisymmetrical and corresponding function $f c {q) 
looks like 

*c(g) = - / sin 2 (7r&(jcos0)d( ? !> = 2-2Jo(27r&(j). (9) 

T J-n 

Although it is not explicitly indicated, functions ^(q) 
depend on the baseline b. Integration over the dimension- 
less frequency q leads to formulae depending on this only 
parameter. Let us denote these integrals as IC(b): 

/>oc 

K(b)= q- 2/3 A(q)9(q)dq. (10) 
Jo 

Integrals fC(b) can be expressed analytically through hyper- 
geometric functions, but exact expressions can not be used 
due to its bad convergence. It is more practical to utilize var- 
ious approximations (see e.g. Martin 1987; Sarazin & Rod- 
dier 1990; Tokovinin 2002). For the c-motion, K(b) is given 
by: 

fC c (b) ^3.367 (1.4023 - b~ 1/3 ), (11) 

Accuracy of this formula is reasonably good because hyper- 
geometric function (entering the precise expression) has a 
good convergence for b > 1 and deviates from 1 by no more 
than 0.005 in domain 1.2 < b < oo. Analogous expression 
for the case of z-tilt can be obtained as well: 

£ c (&)~ 3.367 (1.5024 -b~ 1/s ). (12) 

Formulae (11) and (12) coincide with the equations (7) 
and (8) are given in (Tokovinin 2002) , except that our defini- 
tion of K,(b) differs from the instrument constant K derived 
there in additional factor 6.95 = (0.0229 • 27r) _1 . Also, it 
must be taken into account that K, c {b) = IC t (b) + JCi (b). 

3 WIND EFFECT ON THE DIFFERENTIAL 
IMAGE MOTION 

3.1 Spectral filters for non-zero exposures 

Effect of wind smoothing can be estimated by multiplying of 
temporal power spectrum of the image motion by the spec- 
tral filter corresponding to averaging during the exposure r. 
However, it is easier to use spatial averaging in the direc- 
tion of the wind in rectangular window with length equal to 
the wind shear vt and zero width as it was done by Martin 
(1987). In the used coordinate system, the wind blows at the 
angle 9 to the x axis. Introducing dimensionless wind shear 
u> — vt/D, one can write the spatial averaging filter S(q, (j>) 
as 

Si(q, 4>,lu,9) = sine (u> qcos(4> — 8)), (13) 

where sinc(a;) = sin(7rx)/7ra;. It is clear that the role of this 
filter is to suppress high-frequencies in the spatial spectrum. 

Covariance of image motion between adjacent exposures 
can be expressed through variance of binned data. In order 
to obtain this variance, it is not sufficient just to double 



exposure in formula (13) because in actual measurements 
existing technical limitations lead to a certain pause be- 
tween adjacent exposures. Therefore one should use the fil- 
ter, which corresponds to two separate averaging windows: 

S2(q, <t>,0J,ui,8) = sinc 2 (wflcos(</> — 9)) cos 2 (nuiq cos(</> — 6)), 

(14) 

where the period of exposures T ^ r and corresponding 
dimensionless wind shear Qj — vT/D. It can be seen that in 
the case of T = r expression for £2(0;, o>) becomes identical 
to 5i(2w). 

Multiplication of spectrum F(q, (f>) by the filter S and 
subsequent integration lead to the expression for the vari- 
ance a 2 of the smoothed image motion. Extracting term for 
non-smoothed power oq, one can get 

= F(q,4>)S(q,4>,u,0)dq d<j> 

= 0-0- J I F(fl,0)(l-S(fl,0, W ,0))dfld</>. (15) 

The last integral represents the correction to zero exposure 
5a 2 {uj,9). Hereinafter we will call this value as wind effect 
on variance. 

With help of equation for covariance between adjacent 
positional measurements covi = 2o\ — a 2 , the corresponding 
filter C\ can be obtained: 

Ci(q, <j>,u,u,0) 

= sinc 2 (cj qcos(<f) — 9)) cos(27ro> qcos((j) — 9)). (16) 

Respectively, the wind effect on covariance: 

6c(u,u,0) = jj F(q,4>)(l-Ci(q,<l>,u,u,0))dqd4>. (17) 

The dependence of 8a 2 (ijo,8) and 5c(uj,uj,8) on turbulence 
layer intensity can be removed by normalization by zero ex- 
posure variance cr 2 . 

Parts of the filters Si and C\ depending on angle <f> — 6 
should be included in integrand in expression for function 
*(g) (see section 2.2). Practically, this effect can be esti- 
mated by virtuosic analytics (e.g. by series expansion, Co- 
nan et al. 2000) or numerical integration. The latter is 
more preferable for us since obtained results in any case are 
equally difficult for analysis due to many input parameters. 

3.2 Features of wind smoothing effect 

The dependencies of the normalized wind effects on wind 
direction and wind shear computed with formulae (15) and 
(17) are shown in Fig. 1. Variations of the wind effects with 
wind direction can be clearly seen. The effects are at maxi- 
mum when the direction of image motion coincides with the 
wind direction and at minimum when they are perpendicu- 
lar. The maximum differs from the minimum almost twice. 
On the average, the normalized effects are larger for £-motion 
than for /-motion. The same is valid for the covariance. 

The curves for different values of the wind shear are 
almost similar what indicates that dependencies on 9 and u) 
can be separated. Moreover, when ui < 1 the curves 6c/ 0% 
resemble b~a 2 /o$. Dependence on the wind direction 9 can 
be described by cos 29 with reasonably good accuracy. Note, 
that total variance a 2 of image motion virtually does not 




depend on the wind direction. Its maximal deviation from 
the mean is less than 1.5 per cent. 

Since dependencies of variances erf and of on wind 
direction are 'opposite in phase', their ratio also changes 
markedly. Range of possible ratios erf /erf increases with in- 
crease in wind shear and for u = 1 reaches 1.25 -j- 0.92 of 
initial ratio crf/af. 

The I- and t-motions are independent, in other 
words their correlation (corresponding gradient filter is 
(A/) 2 cos <p sin (f>) equals to zero. However when wind 
smoothing is taking place, the measured I- and t-motions 
become partly correlated. Correlation coefficient changes as 
sin 20, i.e. reaches its maximum when wind blows at angle 
±45° to the baseline. 

The same dependencies as in Fig. 1 but as a functions 
of the wind shear to are presented in Fig. 2. For small argu- 
ment, the curves are parallel to each other what indicates 
their proportionality. Weak dependence of the a 2 on wind 
direction results in a quite small distance between upper 
and lower boundaries. Slope of the curves corresponds to 
quadratic dependence of the effect on u. This feature will 
be considered in details in following section. 



4 SHORT EXPOSURE APPROXIMATION 

4.1 Averaging over all wind directions 

The mean wind effect can be estimated by averaging of the 
expression (15) for the variance a 2 over angle 6. Such aver- 
aging affects only the filter Si and leads to the expression 

° 2 = ff F (<1, <t>) Ti {ojq) dq d<j>, (18) 

where the function Ti(ujq) is a spectral filter of the wind 
shear (Kornilov 2011b). Further averaging over the angle (j> 
results in a new function ^'(q) = 9(q)Ti(u>q). 



The filter of the wind shear has the asymptotic Ti {(jjq) ~ 
1 — TT 2 u) 2 q 2 /6 near the origin which is accurate while co <C l/ir 
because q < 1. Using this approximation, one can calculate 
the wind effect Sa 2 (uj) or correction to zero exposure: 

2 2 poo 

,5a 2 (w) = 2.403AJ D~ 1/3 ^- / q 4/3 A(q)V(q) dq. (19) 

6 Jo 

Integral over dimensionless frequency in this formula 

POO 

C{b) = q 4/3 A(q)y(q)dq, (20) 
Jo 

differs from integral K.(b) in the multiplicand q 4 ^ 3 only and 
also depends on a single parameter b. For the case of c- 
motion, it can be computed analytically and for b = 2.18 
equals to C c {b) — 0.367, varying insignificantly in range be- 
tween 0.374 and 0.366 for 1.2 < b < co. As above, to des- 
ignate the total differential image motion, subscript c are 
used. 

Let us consider the covariance when exposure r = 0. 
Then the filter Si e1 and averaging over all wind directions 
results in 

*'(<?) = *(g) Jo(27rw g ). (21) 

If the wind shear is finite but small, u> -ti 1, an expansion 
of sine 2 near zero can be used and filter Ci(q,<f>,uj,ui,0) can 
be transformed in the following way: 

Ci(q,(j>,u,u),6) = (1- ^ir 2 u> 2 q 2 cos 2 (4>- 6)) 

X cos(27T£2> qcos(<j} — 6>)), (22) 
after averaging this filter over 9: 

*'(«) = *(<?) 

2 2 2 

x (j (27rw g ) - - W & q (Jo(27rwg) - J 2 (2™ q)) J . (23) 
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Figure 2. Normalized wind effect 8cr 2 (ui, 0)/<Jq (top) and 
5c(lo , Co , 0) I CTq (bottom) as a function of wind shear ui and oj. 
For bottom plots ui = 0.8 Co. Upper and lower boundaries for 
total (solid lines), transverse (dashed) and longitudinal (dash- 
dotted) differential motions are presented. The dotted lines depict 
a quadratic dependencies on argument. 



Comparison of this expression with (21) allows to conclude 
that wind effect on covariance depends mostly on the ex- 
posures period T. Even in the case of r = T, the value 8c 
changes only by « 1/6 with respect to the case t — 0. 

Expansion of the Bessel function near zero permits to 
get the formula for 8c: 

8c = +7r 2 a, 2 ) 2.403 AJD- 1/3 C(b). (24) 

From this expression and the formula (19), universal relation 
between the wind effects on variance and covariance can be 
obtained: 

So 2 (u) = 8c(w,w) u2 f^ 2 . (25) 



Table 1. Integrals /C, £, and £* for the b = 2.18. Integrals £* 
for 9 = 7r/4 are equal to ones averaged over wind direction. 



Integral t-motion /-motion c-motion 



g-tilt 



K 




0.8141 


1.3067 


2.1208 


c, 




0.1780 


0.1890 


0.3672 


£*, 


9 = 


0.0489 


0.1402 


0.1891 


£*, 


9 = tt/4 


0.0890 


0.0945 


0.1836 


£*, 


9 = tt/2 


0.1292 


0.0489 


0.1781 








z-tilt 




K 




0.9788 


1.4804 


2.4592 


c, 




0.2206 


0.2314 


0.4518 


£*, 


9 = 


0.0597 


0.1717 


0.2314 


£*, 


9 = vr/4 


0.1103 


0.1157 


0.2259 


£*, 


9 = tt/2 


0.1608 


0.0597 


0.2205 



4.2 Arbitrary wind direction 

For arbitrary wind direction, the computing of the wind ef- 
fects reduces to integration of the expression containing the 
differential filter, the gradient filter and the Si over angle 
<f>. Let us suppose that we can use quadratic expansion of 
I- Sr. 

1 - Si = iyrV q 2 cos 2 (tj>-9). (26) 
Then, the wind effect on variance is calculated by formula 

2 2 

8o 2 {u>,6) = 2.403AJD~ 1/3 ^-r(M), (27) 

where the C* is similar to the integral (20) with functions 
ty*(q) instead ^(q) (see Append. A) and depends not only 
on the b but also on the wind direction 9. 

Values of these integrals C*(b,6) are given in Table 1 
for all components of image motion and for different wind 
directions. Note, that according to (28), the integrals £* 
averaged over all directions are half of the corresponding C. 

Formulae (A4), (A6) and (A7) reflect dependence on the 
wind direction in the form cos 26 as depicted in Fig. 1. These 
formulae allow to deduce expressions for ^*{q,9) averaged 
over all wind directions as well: 

<**M)> = i*(g). (28) 

Therefore, after averaging over 6, formula (27) coincides en- 
tirely with equation (19). 

Similar relation can be obtained for the wind effect on 
covariance by using of quadratic approximation of spectral 
filter 1-Ci: 

2 

1-Ci = +2u 2 ^-K 2 q 2 cos 2 {<t>~9). (29) 

Clearly, the result for Sc(u>, w) differs from expression (27) 
only in factor u; 2 /3 + 2 Cj 2 instead of w 2 /3. It means that for 
arbitrary wind direction, a relation between the wind effects 
on variance and covariance is described by formula (25) as 
well. 
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Figure 3. Accuracy of quadratic correction of variance as a func- 
tion of wind shear ui. Upper and lower boundaries for total (solid 
lines), transverse (dashed) and longitudinal (dash-dotted) differ- 
ential motions are shown. Grey lines correspond to the c-motion 
and coefficient reduced to 0.11 (see the text). Note that for c- 
motion, there are actually two lines for both coefficients but they 
are undistinguished in the picture scale. 

4.3 Applicability of the short exposure 
appr oximat ion 

As it was shown above, the dependence of the wind effects 
on the wind shear demonstrates the quadratic behavior for 
short exposures and period between them. As a matter of 
fact, this is the definition of short exposure regime: the con- 
ditions for which the quadratic approximation is sufficient 
for practical usage. 

In the case of measurement of stellar scintillation with 
Multi- Aperture Scintillation Sensor (Kornilov et al. 2003), 
the effect of wind smoothing can be described by minimal set 
of parameters and a single criterion of applicability of short 
exposure approximation can be defined (Kornilov 2011a). 
For differential image motion, the situation is more complex. 
Fig. 2 indicates that for different components of the image 
motion and wind directions, criteria should be different. 

Acceptable error of the approximation depends on prob- 
lem at hand. As for correction to zero exposure, required 
absolute error of the normalized wind effect coincides with 
desirable relative accuracy e of the result: ASa 2 /^ — e. 

Deviations of the quadratic approximation of &cj 2 ja% 
from their exact functions are shown in Fig. 3. For the l- 
and t-motions, there are a large uncertainty in estimation 
of allowable wind shear. For instance, in order to provide 
accuracy 0.02 of ctq m case °f wind blowing perpendicularly 
to baseline, the wind shear should be less than 0.4 for the t- 
motion and less than 0.8 for the Z-motion. This corresponds 
to wind speed v = 9 and 18ms -1 , respectively. 

Far greater certainty can be reached when the c-motion 
is considered, because it demonstrates weak dependence on 
the wind direction. In this case one can argue that required 
accuracy is provided for w < 0.5 (11ms -1 ) for any wind 
direction. 

The quadratic approximation of the wind effect on co- 
variance Sc/a 2 , works worse, what is expectable from more 
complex behavior of exact dependence and much greater 



coefficient in corresponding formulae. Requirements for ap- 
proximation accuracy will be formulated later in sections 
dealing with practical application of dc/a 2 . 

Quadratic approximation coefficients were computed 
from formula (27) using equation (6) for normalization and 
integrals from Table 1. They ensure correct behavior near 
u) — but do not provide minimal residual in desirable 
range of the wind shear. The range of approximation can 
be expanded by slight reduction of computed coefficient as 
is shown in the figure with the grey curve. For the c-motion, 
error of 5a 2 becomes less than O.OIcto m the range < lu < 1 
which corresponds to maximal wind speed v = 22.5 ms -1 for 
4 ms exposure and 36 ms -1 for 2.5 ms. 



5 EFFECT OF THE WHOLE ATMOSPHERE 
5.1 Total wind smoothing effect 

In previous sections we considered smoothing of the differ- 
ential image motion for infinitely thin turbulent layer with 
intensity A J = C 2 (h) dh and certain wind speed and direc- 
tion. As was mentioned before (see Sect. 2.1), the effect of 
wave propagation was not taken into account, though for 
some cases its impact can be significant 1 . Weak perturba- 
tions approximation assumes independence of phase distor- 
tions generated by different turbulent layers. Because of this 
all second moments can be simply added layer-wise and for 
measurement with zero exposure: 

<7o = 2.403 D~ 1/:i K;(b) [ C 2 (h)dh = 2A03D- 1/3 )CJtot, 

J A 

(30) 

where integration is performed over the whole atmosphere 
and Jtot is the total turbulence intensity. This is the basic 
equation of the differential image motion method of mea- 
surement of the integral turbulence intensity in the atmo- 
sphere. 

For non-zero exposure (or uj 7^ 0) expression is some- 
what more complex because in general case it is necessary 
to take into account dependencies ui(h) and 9{h). Neverthe- 
less, the smoothed variances cr 2 {h) and covariances can be 
integrated over the atmosphere. The same is valid for wind 
effects on the variance Sa 2 and covariance Sc. 

Within the short exposure approximation, total effect 
on the variance (for clarity, dependence on b is omitted) is 
given by: 

8a 2 = 2.403 D- 1/3 ^- f u(h) 2 C 2 (h) C (9) Ah. (31) 

•i J A 

In order to make integrand dependent only on atmospheric 
characteristics, we replace the dimensionless wind shear by 
its definition u> = vt/D and take all device parameters out of 
the integral. Additionally, let us take out value (£*} = C/2 
averaged over all wind directions leaving under the integral 
ratio 2£*(9)/£ which is of the order of unity. 

After these transformations, the integral represents sec- 
ond atmospheric moment of the wind speed distorted by 
angular dependence. Being normalized by total turbulence 
intensity J to t it becomes some effective squared wind speed 

1 http://curl.sai.msu.ru/mass/download/doc/dimm_spccs.pdf 
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(v 2 ). Dividing both parts of expression (31) by Oo, one can 
get 



where 



{v 2 ) = 



Sa 2 



— f 

Jtot J A 



D 2 6 K 



C 2 n (h)v(h) 



(A 



Ah. 



(32) 



(33) 



Similarly, one can obtain the expression for normalized 
wind effect on covariance: 



(34) 



The dependence (32) can be used for estimation of wind 
effect on variance, if vertical profiles of the OT and the wind 
speed or effective squared wind speed are known. On the 
contrary, formula (34) is useful to estimate atmospheric pa- 
rameter {v 2 } from measured 8c value. 

Since sensitivity of DIMM instrument tC(b) increases 
slowly with b and C(b) does not, it follows from equation 
(32) , that DIMM with the larger baseline is less sensitive to 
wind impact. 



0.08 




Figure 4. Deviation of the exact dependence between normal- 
ized wind effects Sir 2 /cr^ and 8c/ a 2 , from (35) as a function of 
wind shear Cj. Designations arc the same as in Fig. 3. Grey lines 
correspond to c- motion with coefficient increased to 0.15 (see the 
text) . 



5.2 Correction of image motion variance 

Recall that in the papers of Martin (1987); Soules et al. 
(1996), the wind correction was estimated on the basis of 
model assumptions about vertical OT profile C 2 (h) and 
wind speed v(h). Estimations of this kind are useful for 
planning or analysis but are not suitable for reduction of 
real data due to low accuracy. 

The usage of the (v 2 ) seems useful for real data re- 
duction. This value is similar to square averaged wind Vi 
which enters the definition of atmospheric coherence time to 
and sometimes being measured by other methods together 
with DIMM measurements. However, without information 
on dominant wind direction these speeds may differ almost 
twice (see Table 1). 

The method of interlaced exposures and the method of 
covariance represent versions of the same method which does 
not require additional data. 

From formulae (32) and (34) it follows as well as for 
single turbulent layer (25) that for short exposure regime, 
the wind effects on variance and covariance are proportional 
to each other and constant of proportionality a depends only 
on parameters of image acquisition process r and T: 

2 T 2 

5a = 5c T 2 +6T 2 = Sca - ( 35 ) 

It allows us to correct the measured variance to zero expo- 
sure using the value 8c also obtained from measurements. 

Constant of the process a = t 2 /(t 2 + 6T 2 ) < 1/7 be- 
cause T ^ r. Both parts of the equation (35) can be divided 
by measured a 2 . In this case procedure of correction looks 
like: 

4 = * 2 (l + affj . (36) 

Recalling definition 8c = ao—covi (17), one can rewrite 
previous expression in the following form: 



The right hand expression includes only measured values, 
where the p is correlation coefficient of successive measured 
centroids. Pay attention that the p — covi/a 2 differs from 
true correlation coefficient covi/a 2 . in normalization. 

Requirements for accuracy of the correlation coefficient 
are quite mild because a < 1/7 (for our device a = 0.0964). 
This is valid not only for measurement errors but also for 
errors arising from quadratic approximation of 8c(Cb)/a 2 . 
Thus, the error 20 per cent in correlation coefficient leads to 
error less than 2 per cent in Co- 
in the method of interlaced exposures, one actually 
measures the variances a 2 and a 2 T affected by wind smooth- 
ing. The short exposure approximation allows to deduce sim- 
ple relation between them: 

2 4 _ 2 1-2 / qo n 

"o = gf T - -a 2r . (38) 

Deviation of the exact dependencies of 8a 2 /a 2 , on Sc/ctq 
from linear functions with theoretical slope are shown in 
Fig. 4. It can be seen that differences are systematically 
positive and grow quite fast. For a 2 /a 2 , error 0.02 is reached 
at w = 0.67. 

For shorter exposure 2.5 ms (instead of 4 ms), the er- 
ror of the approximation significantly reduces (by factor of 
~ 2.5) as well as wind smoothing effect itself. Maximal al- 
lowable value of the wind shear becomes to = 0.9 (16 ms -1 ). 
Because of the constant a becomes 0.04, the impact of 8c 
uncertainty also decreases. 

Slope of the linear dependency predicted by theory, 
does not guarantee the minimal approximation errors in the 
needed range of wind shear. These errors can be significantly 
reduced by increase of slope of linear dependency. For ex- 
ample, the linear approximation of the dependence 8a 2 on 
8c with slope 0.15 provides error less than 0.01 ctq up to 
Cj — 0.8 (14.4ms _1 ). Note, that for our measurements, the 
most probable value Cj = 0.25 and median value is 0.32. 

More complex approximation describing computed de- 
pendence 8a 2 (p) better, does not ensure additivity of the 
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impacts of an individual layers and, in general case, may 
give worse results. Only linear combinations of variance and 
covariance provide independence on OT and wind speed ver- 
tical profile. 



5.3 Wind speed V_ evaluation 

From the expression (33) it follows that (v } differs 
from square averaged wind V2 in additional multiplicand 
2£*(9)/£ which varies by 2-3 times depending on wind di- 
rection in the case of t- and Z-motions. At the same time, for 
the c-motion this factor deviates from 1 by no more than 3 
per cent (see Table f) and angular dependence can be ne- 
glected. Therefore, in this case, one can assume V2 = {v 2 } 1 ^ 2 
exactly. 

The definition of 6c and the expression (37) set the con- 
nection between normalized correction for covariance and 
measured correlation coefficient p: 



what allows to express square averaged wind speed V_ 
through instrumental parameters and measured p: 

^=M^-^T- (4o) 

This dependence gives reasonably good estimation of 
the V2 while quadratic approximation for <5c(w) is valid. For 
greater wind speeds, rise of 5c(Q) slows down and evalu- 
ated V2 becomes underestimated. This effect is illustrated in 
Fig. 5. For our DIMM device substitution of numerical val- 
ues in formula (40) gives V 2 = 13.8 (l-p) 1/2 . Corresponding 
dependence is shown in figure by dashed line. 

Comparison of the exact dependence with approximate 
shows that short exposure approximation underestimates 
wind speed by 13 per cent and 32 per cent for 5ms _1 and 
10ms _1 , respectively. Exposure time has almost no effect 
on appearance of this curve; on the contrary, change in the 
period T leads to proportional change of y axis scale. 

Naturally, impact of fast turbulence layers will be sig- 
nificantly weakened. Therefore, the wind speed derived with 
DIMM can be considered only as lower estimate of the V_. 
To do it more realistic, CCD camera frame period T should 
be reduced twice at least. Then, the underestimate « 30 
per cent will be provided for turbulent layer with 20 ms -1 
wind. Additionally, certain improvement can be achieved by 
increase of DIMM aperture D. Nevertheless, the method is 
able to give quite good estimates even with existent device, 
as will be shown in Sect. 6.3 on real data. 

Considered average wind speed V2 differs slightly from 
V5/3 entering the definition of atmospheric coherence time 
ro. Nevertheless, investigation by Kellerer & Tokovinin 
(2007) has shown that V2 value may be used for estimation 
of to. Recall that method of derivation of ro implemented 
in FADE (Fast Defocusing of stellar image) (Tokovinin, 
Kellerer & Coude du Foresto 2008) and MASS instrument 
(Kornilov 2011a) are based on the V2 either. 




Figure 5. Exact (solid line) and derived by (40) (dashed line) 
dependencies of wind speed V2 on (1 — p) 1 / 2 . 



6 EXPERIMENTAL DATA 

6.1 Measurements data and correction for finite 
exposure 

Experimental proof of the basic relations was made us- 
ing MASS/DIMM data obtained on Mount Shatdjatmaz in 
2007-2010 (Kornilov et al. 2010) for characterisation of OT. 

The program dimm 2 records measurements results in 
output file every 2 s and also its 1-minute averages. These 
results include variances var and covariances of the differen- 
tial image motions and coordinate noise var* in longitudinal 
and transverse directions. Value of coordinate noise is com- 
puted during processing of frames using formula (8) from 
(Kornilov et al. 2007) which gives slightly underestimated 
values for thresholding centroiding method since it does not 
take into account fluctuations of number of pixels exceeding 
a threshold. 

Though in the real conditions, impact of the noise is 
very small (median of ratio var* /var w 0.0008 and in 99.7 
per cent of cases it is less than 0.01), nevertheless, it is 
taken into account in further processing in obvious way: 
a 2 = var — var*. Covariance is not biased by the noise but 
incorrect account for the noise in a 2 can lead to bias in the 
correlation coefficient. 

Distributions of the correlation coefficients in 2007-2009 
(containing « 90 000 1-minute points) have similar charac- 
teristics for all image motions: medians are equal to 0.85, 
only for 4 per cent of events p < 0.6 and for 1 per cent — 
p < 0.44. On the other hand, only for 5 per cent of the mea- 
surements, p > 0.95 and 0.5 per cent > 0.97. Distributions 
of 2-s points demonstrate similar behavior slightly different 
in details. 

Prior to theoretical analysis described in this paper, 
during processing MASS/DIMM data in 2010, we used the 
following empirical law (based on numerical calculations, see 



2 http: / /curl.sai. msu.ru/mass/download/doc/infrasoft_eng.pdf 
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Figure 6. Distribution of normalized corrections <5o- 2 /<7g com- 
puted using formula (37) with coefficients a = 0.096 (solid line) 
and a = 0.13 (dashed line). 



Sect. 3.1): 



(To = ^ (1 + 0.15(1 -pj). 



(41) 



To define the coefficient 0.15, the differential distribution of 
the p as weights was used. Dependence (41) corresponds to 
the formula (37) with a = 0.13. 

In Fig. 6, the distributions of the corrections da 2 /oq ob- 
tained from the p using formula (37) with theoretical values 
a = 0.096 and a = 0.13 are shown. These distributions are 
quite close. Median of required correction is about 0.02. For 
90 per cent of measurements, the correction is less than 0.04. 
It is possible to estimate that in marginal cases p < 0.5, the 
correction reaches 0.08. Note that these numbers are quite 
expectable taking into account that site under consideration 
has small median wind speed (2.3 ms -1 ). 



6.2 Dependence on wind direction 

Analysis of angles between wind direction and baseline for 
our measurements shows that their distribution is close to 
uniform. There is no any preferable angles, variations does 
not exceed ±20 per cent. For the detection of angular depen- 
dencies in the I- and the Amotion, the following combination 
of the wind effects are used: 



e 



6c t (0) -5ci(6) pi-pt 



Sc c (9) 1 — p c 

Using results of Sect. 5.1, it can be shown that 

{cue) -cm) 



6 : 



vm) 



(42) 



(43) 



where averaging is done over the whole atmosphere with 
weight v 2 C 2 . Assuming that for majority of the measure- 
ments the ground layer is dominant, one can use surface 
wind direction as estimate of mean wind direction. 

In Fig. 7, computed dependencies of O ratio on wind di- 
rection for different values of w (they affect O only slightly) 
are shown. Medians, 1st and 3rd quartiles (in groups of 500 
points) of 6 derived from measured p are also presented. To 




180 



Figure 7. Median © in groups of 500 points as a function of 
angle between wind direction and baseline (circles). The 1st and 
3rd quartiles are indicated by dashed lines. Thin lines stand for 
expected behavior of for different values of ui 



minimize the impact of measurements errors and uncertain- 
ties of the wind direction, data points were selected when 
surface wind was vgl ^ 4ms -1 . 

Observed dependence looks similar to expected though 
its amplitude is lower. Difference in amplitude can be ex- 
plained by the fact that wind direction can be different in 
different layers of atmosphere. Also amplitude could be re- 
duced by small negative bias in p c . We did not find depen- 
dence of the correlation coefficient of the c-motion on the 
wind direction. 



6.3 Comparison of wind speeds estimated from 
MASS and DIMM data 

For verification of the estimate the effective wind, we have 
performed comparison of V2 estimates from data obtained 
simultaneously with DIMM and MASS. In order to do this, 



the average wind speed Vf ree = 1-59 -10 J t 



-3/5 

free 



t„ 1 in free 

atmosphere was computed using to and free atmosphere in- 
tensity Jf rco measured with MASS (Kornilov 2011a). Then, 
the wind speed was corrected with ground layer intensity 
Jgl and surface wind speed «gl from meteodata with help 
of expression V^ass = (Ffr«Aee + ^gi^gl) / 'Jtot- 

The frame rate of our DIMM camera reaches 200 frames 
per second. The program dimm output data format allows 
us to compute p c and then V2 using formula (40). 

Results of comparison of 1-minute V2 estimates and 
wind speeds Umass is shown in Fig. 8, in the form of me- 
dian and 1st and 3rd quartiles (in groups of 500 points). The 
quartiles indicate that the dependence is quite narrow, mean 
spread is ±1 ms _1 . Median curve is similar to expected ex- 
act dependence depicted by thin solid line but for higher 
wind speed it lies somewhat lower (this reflect the loss of 
impact of fast turbulence layers). For clarity, left vertical 
axis presents the atmospheric time computed assuming me- 
dian fio = 0.93 arcsec. 

Median values of the wind speeds Vmass = 5.84 ms -1 
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Figure 8. Comparison of V2 estimated from DIMM and MASS 
data: dots stand for median wind speed in groups containing 500 
measurements, dashed lines stand for 1st and 3rd quartiles of 
distributions in these groups. Solid line stands for expected de- 
pendence considering errors of short exposure approximation. 



and V2 = 5.67 ms" 1 are very close, though there are reasons 
for believing that surface wind speed vgl is underestimated 
due to anemometer was located at altitude 5 m. Nevertheless 
in median region wind speeds coincide quite well, in spite of 
absence of any empirical calibrations. 



7 DISCUSSION AND CONCLUSIONS 

Applicability of quadratic approximation of the wind effects 
for the measurement of the differential image motion de- 
pends on several instrumental parameters. Exposure time 
and period between exposures are most important of them. 

Data analysis indicated that for two brightest stars used 
in our measurements, 4 ms exposure lead to saturation in 
centres of images. Because of this we reduced exposure to 
2.5 ms starting from December, 2009. Median noise impact 
increased insignificantly to w 0.0011 and for the faintest 
target stars with magnitude ~ 3 m became « 0.004. Such 
negligible increase of the noise means that exposure can be 
reduced further. 

The 2.5 ms exposure requires correction less by factor 
of approximately 2.5 (a = 0.04). For our measurements it 
means that in 90 per cent of time correction is less than 
0.015 and only in rare cases ~ 0.033. 

Reduction of the exposures period T considerably ex- 
pands capabilities of the DIMM. For the most CCD, maxi- 
mal frame rate is defined by exposure and height of region 
of interest on detector. Usually, the height of the window is 
set to 60-80 CCD rows. It is conditioned by amplitude of 
telescope tube vibrations occurring at wind speed greater 
than ~ 7ms" 1 . More robust mount will allow to reduce the 
height of the window to 20-30 rows and increase frame rate 
up to 400 frames per second. 

Such acquisition rate makes it possible to estimate mean 



squared wind V2 with reasonably good accuracy up to speed 
~ 10 ms" 1 , what corresponds to to = 3 ms for site with 
median seeing /3o = 1 arcsec. Lower values of to will be esti- 
mated with larger bias but these situations are less interest- 
ing for high-technology astronomical observations. Anyway, 
the possibility to estimate the atmospheric coherence time 
using DIMM is promising because DIMM is most common 
site testing instrument. 

Let us briefly summarize main conclusions from anal- 
ysis of the differential image motion in the short exposure 
approximation given in the paper. 



(i) Transverse image motion is more subject to the wind 
smoothing. The effect is at its maximum at wind blowing 
across the device baseline. On the contrary, for i-motion, 
the effect is maximal at wind blowing along the baseline. 

(ii) Total differential image motion o 2 — of + of has a 
simpler description. Its main advantage is a weak depen- 
dence (less than ±0.02) on the wind direction. 

(iii) Effect of the wind smoothing for short exposures is 
proportional to exposure squared. It is possible to signifi- 
cantly diminish effect without loss of accuracy by using a 
shorter exposures. 

(iv) Effect of the wind smoothing for short exposures de- 
pends on device aperture diameter as oc -D" 2 , what leads 
to significant increase in wind effect for DIMMs with small 
apertures. 

(v) Correction of measured variance of image motion can 
be performed in simple way using measured correlation co- 
efficients between adjacent frames. 

(vi) Value of the correlation coefficient for c-motion does 
not practically depends on wind directions and can be used 
for the estimation of the averaged wind V2 and the atmo- 
spheric coherence time for the whole atmosphere. 



Taking into account the fact that power of the c-motion 
a 2 can be derived with greater accuracy than for of and 
of individually, we recommend to use it and corresponding 
relation (11) or (12) for calculation of the integral turbulence 
intensity. Widely spread estimation of the seeing as average 
between /3o obtained from /- and i-motions in conditions of 
significant wind smoothing is biased and non-optimal. 

Expressions for the wind effect Sc obtained in this work 
can be adopted to another pattern of readout from DIMM 
camera. For example, with full-frame CCD it is possible to 
obtain 2 (or more) successive exposures without pause be- 
tween them and then to readout all subframes together dur- 
ing some time. 

In the present work we did not consider influence of in- 
ner and outer scale of turbulence and possible departures 
from Kolmogorov law. We also did not take into account 
the effect of propagation of light wave which can be sig- 
nificant for high-altitude turbulence. Effect of propagation 
suppresses high frequencies of spatial spectrum of the im- 
age motions and, as a consequence, an influence of the wind 
smoothing become smaller. This reflects in decrease of inte- 
grals K, and C, and the latter changes far more significantly. 
Amount of effect is about 20 per cent for the turbulence in 
the tropopause. Nevertheless, the quadratic approximation 
(32, 34) is still valid and can be used. 
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APPENDIX A: * FUNCTIONS 

Averaging of components of differential motion spectrum de- 
pending on angle cj> over this angle results in the expressions 



2 f+' r 

*;(?) = — sin (Tvbqcoscf)) cos <fid(f> 



1 - Jo(2nbq) + J 2 (2TTbq), (Al) 



and 



2 r +1T 

^t{q) = — I sin (irbq cos </>) sin (f>dcj> 

= 1 - Jo (271-6(7) - J 2 (2-Kbq). (A2) 



Here the integrals J7o and J\ from Append. B were used. 

For c- motion er 2 = of + of , result of the averaging is 
equal to sum *;(</) + ^t(q) or 



*c(«) 



sin (irbq cos <j>) d<j> = 2 - 2 Jo (2nbq) . (A3) 



In presence of wind smoothing, a similar functions 
ty*(q,8) can be obtained. In order to do this one should 
multiply integrands by cos (<p — 9) what follows from (26) 
and (29). This factor can be represented as cos 2 (0 — 9) = 
cos 2 9 — cos 29 sin 2 (f>. The cross term with cos <j> sin <fi is omit- 
ted because it gives zero after integration. Using formulae 
from Append. B it is possible to write: 

K(q) = cos 2 9 * ; (q) + sin 2 9 * t (q) 

= 1 - Jo {2nbq) + cos 29 J 2 (2-Kbq) . ( A4) 

For wind blowing across the DIMM base, function **(?) 
coincides with ^i(q) and for wind along, with vE , t(q). For 
other wind directions ^c(q) lies between these boundary 
functions. 

For other components we will need the value of the in- 
tegral: 



^ 2 (q) = — / sin {-nbq cos <f>) sin 

7T 



3 6J 2 (2Trbq) 
4 



(2nbq) 2 
(A5) 

Using relation **(g) = &*{<}) ~ **(?), the following expres- 
sions for I- and t-motions can be obtained: 

*«*(«) = ^MQ) + cos2 (5*1(9) - Ml) + *»(«))■ (A6) 
***(<?) = \*t{q) + cos 20 (±* t fo) - * 2 (g)). (A7) 



APPENDIX B: BASIC INTEGRALS 

The expressions from Append. A reduce to tabular integral 
3.715.21 (Gradstein & Rygik 1963): 

/•7T 

2„(z) = / cos(0cosa:) sin 2 " xdx 
Jo 

= ^r(n + l/2)(2/z) n J n (z). (Bl) 

In calculations we used integrals with n = 0, 1,2 in limits 
[— 7T, +7r]: Xo{z) — 2iv Jo(z), Ii{z) = 27T Ji(z)/z and Ti(z) = 
67T J2(z)/z 2 in following form (as a rule): 



J, 



/ + 7T 
sin 2 {z cos x) 
-TV 

J — 7T 



sin n xdx 



cos(2zcos:r)) sin " xdx. (B2) 



In particular, J (z) = tt(1 - J (2z)), Ji(z) = n(l/2 - 
J 1 (2z)/2z) and J 2 (z) = tt(3/8 - 3 J 2 (2z)/4 2 2 ). Integrals 
with cos 2n (a;) term obviously reduce to ^7-integrals. Inte- 
grals with sin to the odd power in mentioned limits are equal 
to zero. 



